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We study various statistical properties of real roots of three different classes of random polynomials 
which recently attracted a vivid interest in the context of probability theory and quantum chaos. 
We first focus on gap probabilities on the real axis, i. e. the probability that these polynomials have 
no real root in a given interval. For generalized Kac polynomials, indexed by an integer d, of large 
degree n, one finds that the probability of no real root in the interval [0, 1] decays as a power law 
n~ s ' d ' where 6(d) > is the persistence exponent of the diffusion equation with random initial 
conditions in spatial dimension d. For n ^> 1 even, the probability that they have no real root 
on the full real axis decays like n~ 2 ^ e( - 2 ^ +e< - d ^ . For Weyl polynomials and Binomial polynomials, 
this probability decays respectively like exp (— 2# 00 y / n) and exp (— ■K0oo\/n) where 9oo is such that 
9(d) = 2 _3/ ' 2 # 0O \/d in large dimension d. We also show that the probability that such polynomials 
have exactly k roots on a given interval [a, b] has a scaling form given by exp (—N a b (p(k/N a b)) where 
N a t is the mean number of real roots in [o, 6] and <p(x) a universal scaling function. We develop 
a simple Mean Field (MF) theory reproducing qualitatively these scaling behaviors, and improve 
systematically this MF approach using the method of persistence with partial survival, which in 
some cases yields exact results. Finally, we show that the probability density function of the largest 
absolute value of the real roots has a universal algebraic tail with exponent —2. These analytical 
results are confirmed by detailed numerical computations. Some of these results were announced in 
a recent letter [G. Schehr and S. N. Majumdar, Phys. Rev. Lett. 99, 060603 (2007)]. 

PACS numbers: 02.50.-r, 05.40.-a,05.70.Ln, 82.40.Bj 

I. INTRODUCTION 

Despite several decades of research, understanding the zero crossing properties of non-Markov stochastic processes 
remains a challenging issue. Among them, the persistence probability pq (t) received a particular attention, especially in 
the context of many-body non-equilibrium statistical physics, both analytically [l[ as well as experimentally @, S 0, [f| . 
The persistence po (t) for a time dependent stochastic process with zero mean is defined as the probability that it has 
not changed sign up to time t. In various physical situations, po(t) has a power law tail po[t) ~ t~ e where 9 turns out 
to be a non-trivial exponent whenever the stochastic process under study has a non Markovian dynamics. One such 
example is the diffusion, or heat equation in space dimension d where a scalar field 0(x, t) evolves according to the 
deterministic equation 

9 t 0(x,<) = V 2 0(x,t), (1) 

with random initial conditions <fi(x,t — 0) = tp(x) where il>(x) is a Gaussian random field of zero mean with delta 
correlations [?/;(x)V>( x ')]mi = 8 d (x — x'). We use the notation [...];„; to denote an average over the initial condition. 
For a system of linear size L, the persistence Po(t, L) is the probability that 0(x, t), at some fixed point x in space, 
does not change sign up to time t. The initial condition being (statistically) invariant under translation in space, this 
probability does not depend on the position x. In the scaling limit f > 1, L > 1 keeping the ratio t/L 2 fixed, it was 
found in Ref. [6| that po (t, L) takes the scaling form 



Po 



(t,L) cx L- 26 ^h(L 2 /t) , (2) 



where h(u) ~ c st , a constant independent of L and t, for u <C 1 and h(u) cx u 6 ^ for u ^ 1 where 8(d) is a 
d-dependent exponent. This implies that in the L — > oo limit, po(t) = po(t,L — > oo) ~ for large t. It 

was shown in Ref. @ that the probability Vq(T) that a Gaussian stationary process (GSP) with zero mean and 



2 



correlations [cosh(T/2)]~ d / 2 decays for large T as Vq(T) ~ exp [— 9{d)T] where 9(d) is the same as the persistence 
exponent in diffusion equation. This exponent 9(d) was measured in numerical simulations @,0|j yielding for instance 
#sim(l) = 0.12050(5), #sim(2) = 0.1875(1). The case of dimension d = 1 is particularly interesting because 9(1) was 
determined experimentally using NMR techniques to measure the magnetization of spin polarized Xe gas [B[ , yielding 
#exp(l) = 0.12 in good agreement with numerical simulations. In the limit of large dimension d, which will be of 
interest in the following, one can show that 9(d) = 2~ 3 / 2 # 00 v / d where 9 X is the decay constant associated with the 
no zero crossing probability of the GSP with Gaussian correlations exp (— T 2 /2), which was studied in the past by 
engineers, in particular in the context of fading of long- wave radio signals (see for instance Ref. Q). 

A seemingly unrelated topic concerns the study of random algebraic equations which, since the first work by 
Bloch and Polya 0] in the 30's, has now a l ong story [l(| EH- Recently it has attracted a renewed interest in 
the context of probability and number theory [k| as well as in the field of quantum chaos [l3j]. In a recent letter 
(l4| . we have established a close connection between zero crossing properties of the diffusion equation with random 
initial conditions U) and the real roots of real random polynomials (i.e. polynomials with real random coefficients). 
In Ref. [3], we focused on a class of real random polynomials K n (x) of degree n, the so called generalized Kac 
polynomials, indexed by an integer d 



d-2 



K n (x) = a Q + ^ai i 4 x l . (3) 



i=i 



Here, and in the following, dj's are independent real Gaussian random variables of zero mean and unit variance 
(aidj) = Sij where we use the notation (...) to denote an ave rag e over the random coefficients a^. In the case of d = 2, 
these polynomials reduce to the standard Kac polynomials [15j |. which have been extensively studied in the past (see 
for instance Ref. [12] for a recent review). In that case, we will see below that the statistics of real roots of K n (x) is 
identical in the 4 sub-intervals [— oo, —1], [—1,0], [0, 1] and [1, +oo]. Instead, for d ^ 2, which was studied in Ref. [lq |. 
the statistical behavior of real roots of K n (x) depend on d in the inner intervals, while it is identical to the case 
d = 2 m the outer ones. Focusing on the interval [0, 1], we asked the question : what is the probability Po([0, x], n), 
< x < 1, that K n (x) has no real root in the interval [0,x] ? Such probabilities were often studied in the context of 
random matrices, where they are known as gap probabilities (l7| and in a recent work [l8| . Dembo et al. showed that, 
for random polynomials K n (x) with d — 2, Po([0, l],n) oc n - ^ 2 ) where the exponent ((2) = 0.190(8) was computed 
numerically. In Ref. [l4|| . by mapping these two random processes (fT]) and ([3]) onto the same GSP, we showed that in 
the limit 1 - i « 1 « n keeping n(l — x) fixed, one has (similarly to Eq. |(3J)) 

Pq([0, x], n) oc rC e{d ^h-(n(l -x)) , (4) 

with h~(y) — ► 1 for y -C 1 and h~(y) ~ y 8 ^ for y 1, yielding in particular Po([0, i\,n) oc nT^^ thus identifying 
C(d) = 9(d). We then extended our study to the probability Pk([0, 1], n) that generalized Kac polynomials have 
exactly k real roots in [0, 1] and we showed that it has an unusual scaling form (for large k, large n, but keeping the 
ratio k / log n fixed) 

./ k 

M 



P fe ([0,l],n)ocn ^Vi°s«; , (5) 

where <p(y) is a large deviation function, with <fi(0) = 9(d). In both cases, our numerical analysis suggested that h~ (y) 
and (p(y) are universal in the sense that they are independent of the distribution of provided (a 2 ) is finite. The 
purpose of the present paper is twofold : (i) we will give detailed derivations of the results announced in Ref. [l4| 
together with some new results, like the distribution of the largest real root, for generalized Kac polynomials K n (x) ([3]), 
(ii) we extend these results (j4] [5]) to two other classes of random polynomials which were recently considered in the 
literature. First we will study Weyl polynomials W n (x) defined as 

W n (x) = J2 a i^- (6) 

i=0 * ! ' 



Recently, the distribution of complex zeros of Weyl polynomials with complex coefficients were observed experimentally 
in a degenerate rotating quasi- ideal atomic Bose gas [19(. Here we will focus instead on the real roots of such 
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polynomials with real coefficients. Besides, we will consider binomial polynomials B n (x) denned as 



i=0 



As is pointed out by Edelman and Kostlan [12J, "this particular random polynomial is probably the more natural 
definition of a random polynomial" . In the literature, they are sometimes called SO (2) random polynomials because 
their m-point joint probability distribution of zeros is 5*0(2) invariant for all m [20|. We will show below that the 
gap probabilities for these classes of random polynomials ([HI [7J are closely related to the persistence probability for 
the diffusion equation in the limit of large dimension. Our main results, together with the layout of the paper, are 
summarized below. 



1. In section II, we briefly recall the main properties of the persistence probability, po(t,L), for the diffusion 
equation with random initial conditions. In section II- A, we recall the finite size scaling for po(t, L) in dimension 
d whereas in section II-B, we focus on the limit d — > oo. 

2. Section III is devoted to real random polynomials, where our main results are presented. In section III-A, 
we present a detailed study of the density of real roots for these three classes of polynomials, which turns 
out to behave quite differently in all the the three cases under investigations ([3J [5J UJ ■ In section III-B, we 
will turn to the analysis of gap probabilities, which we will first analyse from the point of view of two-point 
correlations. Next, we will present a mean field approach, or Poissonian approximation, which neglects the 
correlations between the real roots of these polynomials, to compute the gap probabilities. We will further show 
how this mean field approximation can be systematically improved using the persistence probability with partial 
survival [2l|, which in some cases even yields exact results. In particular we show that the probability qo(n) 
that these polynomials have no root on the full real axis is given by 

q (n) ~ n -md)+e(2)) iovK n (x), 
q (n) ~ exp (-26*00 y/n) for W n (x) , 

q (n) ~ exp (-nOooy/ri) for B n (x) . (8) 

In section III-C, we will then generalize this study to the probability that these polynomials have exactly k real 
roots on a given real interval. Extending the results obtained in Ref. [3 for K n (x) like in Eq. ©, to Weyl and 
Binomial polynomials, we will show that the probability qk{n) that W n (x) and B n (x) have exactly k roots on 
the full real axis has a scaling form (for large k, large n, but keeping the ratio k/y/n fixed) given by 

q k (n) ~ exp [-y/nip (k/y/n)] , (9) 

where <p(y) is a large deviation function, which depends on the polynomials under consideration W n (x) or B n (x). 
We will also show that these scaling forms in Eq. ([5l [9]) can be qualitatively described by the aforementioned 
mean field approximation. To end up, we study in section III-D the probability density (p.d.f.) p m ax(a;) of the 
largest absolute value of the real roots and obtain the exact asymptotic result 

Pmax(x) OC \ , X > 1 , (10) 
X 

for all the three classes of random polynomials under investigation. All our analytical results will be verified by 
numerical computations and some details of the analytical computations involved in this section have been left 
in Appendices A,B, C, D and E. 

3. Finally section IV contains our conclusions and perspectives. 
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II. A BRIEF OVERVIEW ON PERSISTENCE FOR DIFFUSION EQUATION IN DIMENSION d 

A. Persistence exponent 8(d) and finite size scaling 

We consider a scalar field 0(x, t) in a d-dimensional space which evolves in time under the diffusion equation (fTj). 
For a system of linear size L, the solution of the diffusion equation in the bulk of the system is 

0(x,t)=/ <fy0(x-y,i)V(y) , S(x) = (4^)~ d/2 exp(-x 2 /4t), (11) 

J|y|<L 

where ^>(x) = </>(x, 0) is the initial uncorrelated Gaussian field. Since Eq. pT|) is linear, 0(x, f) is a Gaussian variable 
for all time t > 0. Therefore its zero crossing properties arc completely determined by the two time correlator 
[0(x, t)<p(x, i')]ini- To study the persistence probability po(t,L) it is customary to study the normalized process 
X(t) = 0(x, i)/[0(x, t) 2 ]^; 2 [H- Its autocorrelation function a(t,t') = [X(t)X(t')]i n i is computed straightforwardly 
from the solution in Eq. (fTTj) . One obtains a(t, £') = a(t, £') with i = t/L 2 , t' = t 1 / L 2 and 



ill 1 



d/4 

t, t' < 1 



a(t,t') = { \WWJ ' ^ (12) 
1 , t, t' > 1 . 

We first focus on the time regime t,t' <C 1. In terms of logarithmic time variable T — \ogt, X(T) is a GSP with 
correlator 

a(T, T') = a(T - T') = [cosh(|T - T'\/2)]- d / 2 , (13) 

which decays exponentially for large \T — T'\. Thus the persistence probability p (t,L), for t L 2 , reduces to the 
computation of the probability Vq(T) of no zero crossing of X(T) in the interval [0, T]. It is well known [22j that if 
a(T) < 1/T at large T then Vo(T) ~ exp[— 9T] for large T where the decay constant 9 depends on the full stationary 
correlator a(T). Reverting back to the original time t — e T , one finds po(t,L) ~ t~ e( * d \ for t <C L 2 . In the opposite 
limit t ^> L 2 , one has po(t, L) — > Al, a constant which depends on L. These two limiting behaviors of po(t, L) can be 
combined into a single finite size scaling form as in Eq. 

Despite many efforts, there exists no exact result for 6(d). However various approximation methods have been de- 
veloped to estimate it. One of the most powerful is the so called Independent Interval Approximation (IIA) [23[, 
which assumes the statistical independence of the intervals between successive zeros of (f>(x,t). This gives e.g. 
#iia(1) = 0.1203..., #ha(2) = 0.1862... 6] in remarkable agreement with numerical simulations. A more systematic 
approach is via persistence with partial survival [2lJ, which we will use below (see section III-B). An alternative 
systematic approach consists in performing a small d expansion [24[ yielding 9(d) — d/4 — 0. 12065. ..d z l 2 + which 
would certainly require higher order terms to make it numerically competitive. Yet another systematic approach is a 
series expansion introduced in the context of "discrete time persistence", yielding results for 9(d) which are in very 
good agreement with numerical simulations (25| . 



B. Persistence in the limit of large dimension d 



As we will see later, some statistical properties of the real roots of the polynomials W n (x) © and B n (x) (|7|) turn 
out to be related to the statistics of zero crossings of the diffusion equation in the limit of large dimension d. To study 
the persistence probability in that limit one performs a rescaling of the T variable in Eq. (fT5|) . T = 2 3 / 2 T '/yd such 
that 



a(T — T') = a 2 3 / 2 



f-T' 
Vd 



exp 



--(f-T 
2 V 



'\2 



d> 1 



(14) 



Therefore in the limit of large dimension d, one has 9(d) — 2~ 3 l 2 9 oa \fd where 9^ is the decay constant associated 
with the no zero crossing probability of the GSP with correlator exp [— \(T — T') 2 ]. Even in that limit, there is no 
exact result for 9^. However, it can be approximately estimated using IIA [gsj, yielding (^ha = 0.411497... 0| in 
very good agreement with numerical simulations #oo,sim = 0.417(3) 0]. 



■5 



III. RANDOM POLYNOMIALS 

We now focus on statistical properties of the real roots of random polynomials, extending our previous study 
presented in Ref. [3]. Being Gaussian processes, the statistical properties of these polynomials are determined by 
the 2-point correlators C n (x,y), given by 

C n (x,y) = {K n (x)K n (y)) = \ + 2_^i 2 (xy) 1 for Kac polynomials , (15) 

i=l 

CJx, y) = (W n (x)W n (y)) = V ^ for Weyl polynomials , (16) 

C n (x, y) = (B n (x)B n (y)) = ^ ( , j (xy)' 1 = (1 + xy) n for Binomial polynomials . (17) 



i=0 



For simplicity, we chose the same notation C n (x, y) for the three classes of polynomials under study, and we will do 
so for other quantities. In the following, these three polynomials will be treated separately so this should not induce 
any confusion. For later purposes it is convenient to introduce the normalized correlator C n (x, y) with 

C n (x,y) = —= . (18) 

^/C n (x,x)C n (y,y) 

Notice that C n (x,y) — C n (l/x, 1/y) for Kac polynomials K n (x) with d = 2 and for Binomial polynomials B n (x). 

A. Density and mean number of real roots 

Let us denote Ai, A2, — , A p the p real roots (if any) of one of these random polynomials in Eq. ([3j[6j[7]). The mean 
density of real roots p n (x) is given by 

p 

p n (x) = y^(£(x - A^) = (\K' n (x)\5(K n (x))) , for Kac polynomials , (19) 

i=l 

and similarly for Weyl polynomials W n (x) and Binomial polynomials B n (x). Under this form (|19|) . one observes that 
the computation of the mean density involves the joint distribution of the polynomial K n (x) and its derivative K' n (x) 
which is simply a bivariate Gaussian distribution. Thus computing p n (x) involves a double integration of a bivariate 
Gaussian distribution. This can be easily performed to obtain the following well known result 

PnM - ^EMh±<mEE^ , Cn{x) = Cni x, x) . (20) 

Z7TC n yX J 



This formula (|20|) can be written in a very compact way [12I ] : 



P 



„{x) = -\/d u d v \0gCn(u,v) 



(21) 



u—v—x 



For these different polynomials in Eq. (J3j [6l [Tj) , we will be interested in the number of real roots on a given interval 
[a, 6], which we will denote N n [a,b]. Being a random variable, we will focus on its moments (N^[a, b]), with k G N. 
In particular, one has from the definition of p n (x) in Eq. (|19[) 

(N n [a,b]) = [ pn(x)dx , (22) 



and higher cumulants will be considered below. 
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1. Generalized Kac polynomials 

One remarkable property of the generalized Kac polynomials K n (x) is that, in the large n limit, the roots in the 
complex plane tend to accumulate close to the unit circle centered at the origin. In the left panel of Fig. [TJ we show 
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FIG. 1: Left : Mean density of real roots p n (x) for Kac polynomials K n (x) © and d = 2 as a function of x for different values 
of n = 10, 50 (dotted lines). The solid line is the analytic expression poo{x) in Eq. (|25|l for \x\ < 1 and in Eq. I|27p for \x\ > 1. 
Right : Plot of n~ 1 p n (x) as a function of n(l — x) for n — 500, 1000 (and thus x close to ±1). The dotted line is the function 
p K (v) m Eq. Q31|). There is no fitting parameter. 

a plot of the density of real roots p n (x) for d = 2 computed from Eq. (|2"Tj) for different values of n — 10 and 50. In 
the large n limit, one clearly sees that the real roots of such polynomials are concentrated around x = ±1, where the 
density is diverging. This can be seen by computing p n (±l) from Eq. (|21|) 



1 ((1 + Hijh 1 - d/2))H(n, -d/2 - 1) - H(n, -d/2) 2 ) 2 2n 
Pn{±1} ~ n l + W(n,l-d/2) a ir(d + 2) 



(23) 



where H(n, r) — Ylk=i xs a generalized harmonic number [27j • To obtain the asymptotic behavior in the large n 
limit of the above equation (|23|) we used H(n, r) oc / (1 — r), for large n and r < 1. 

Away from these singularities, p n (x) has a good limit when n — > oo. However, one has to treat separately the cases 
\x\ < 1 and \x\ > 1. For |x| < 1, the calculation is straightforward because C n (x,y) in Eq. (I15[) has a good limit 
n — > oo when x, y < 1. This yields 



Poo (a;) 



Li-i_ d / 2 (^ 2 )(l + Lii_ d/2 (x 2 )) - Li 2 



d/2 



(x 2 ) 



Ixl < 1 



^|x|(l+Li 1 _ d/2 (x 2 )) 

where Li„(z) = X^i z V*" i s the polylogarithm function (27}- In particular, one has Poo(0) 



(24) 

l/7r for all d, and 



5 (x) - (d/2)2(27r((l - x)))- 1 for x 
1 



r 



7r(l 



d = 2 



For instance, one has for |x| < 1 
1 



7T(1 



2x 6 - 4x 4 + 2x 2 + 1 
2x 6 + 3x 4 - 2x 2 + 1 



d = 4 



(25) 



For |x| > 1, the analysis is different because the correlator C n (x,y) in Eq. (fT5")) does not converge any more in the 
limit n — > oo when x, y > 1. Instead, one has in that case (see also Ref. |l6|) 



rf-2 

C7„(x, y) = 1 + V i— (xyY oc [ -^>— 



x,y > 1 



(26) 
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This leads to the expression for the density Poo(x) for |x| > 1 : 

Poo(x) = —-^ rr , (27) 

ir(x z — 1) 

independently of d. To understand better the divergence of p n (x) around x = ±1 (|23| when n» 1, we focus on Pn(x) 
around x = 1. In the limit n> 1 and 1 — a; -C 1 keeping y = n(l — x) fixed, one shows in Appendix [Al ( see also [28j) 
that 



Ki n \\ k i \ 1 I d/2+i(y) ( h/2(y) , ,,, • 



J m(2/) = / dx a;™ exp (-2yx) . 
Jo 



(29) 



One has p K (0) = ^ y , recovering the large n behavior in Eq. and its asymptotic behaviors are given 

by (see Appendix lA~f 

-HfV 3 ' ^ + °° (.-so) 



27T|tf| 



y -oo 



1/1 1 \ 1/2 

P K (2/) = ^- - in d = 2. (31) 



For instance, one has 



27r \ y 2 sinh y 

In the right panel of Fig. [1] we show a plot of n~ 1 p n (x) 1 where p n (x) is given in Eq. (|A1[) , as a function of n(l — x) 
for d = 2 and different large values of n = 500, 1000 together with the asymptotic results in Eq. (f3"Tj) : we find a very 
good agreement with these analytic predictions (|28l l3Tj) . 

Having computed the mean density of real roots, we now focus on (N n ([a,b\)}. On the interval [0,1] the main 
contribution to the mean number of real roots comes, for large n, from the vicinity of x — 1. Therefore, to compute 
(N n [0, 1]) to leading order in n, one uses the scaling form for the density in Eq. (|28|) . valid close to x = 1, and the 
asymptotic behavior in Eq. (|30|) to obtain for n ^> 1 



(N n [0,l]) = {N n [-l,0]) = J p K {y)dy + 0{l) = — S j-\ogn + 0{l), (32) 

where the corrections of order 0(1) receive contributions from the whole interval [0, 1] (not only from the vicinity of 
x = ±l). 

Similarly, one gets from Eq. (gSJ) and Eq. f5D]l: 

C n 1 
(JV n [-oo,-l]) = <JV n [l,+oo])= / p K {-y)dy + 0(1) = — logn + 0(1) , (33) 

Jo 27r 
which is independent of d fl6l ] . From Eqs (|32l I33|) we compute the total number of roots on the real axis : 



(JV n ([-oo,+oo])) = -[l + i- 2 ) ^gn + 0(l) , (34) 

thus recovering, in a way similar to the one used in Ref. [12] for d = 2, the result of |l6[. Notice that for d = 2 
the higher order terms of the large n expansion in this formula (|34[) have been obtained by various authors (see for 
instance Ref. (l2l, [29j ) , although, to our knowledge, they have not been computed for d^2. 
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In view of future purposes, we also compute (N n [0,x}} in the asymptotic limit where n ^> 1, and < 1 — x <C 1 
with n(l — x) kept fixed : 

(JV n [0,a:])=(JV B [0,l])-»»-(n(l-a!)) , (35) 

V-(V)= I du P K {u), (36) 



such that r]- (0) = and with the asymptotic behavior obtained from Eq. (|30|) 

v-(v) ~ ir\ o lo sy i y-^oo. (37) 

Z7T V A 

Similarly, we compute (JV n [x, oo]) when x > 1 and obtain for n ^> 1 and 0<i-1<1 keeping y = n(x — 1) fixed 

(N n [x,oo\) = (N n [l,oc})-ri + (n(x-l)) , (38) 

V+ (y) = f J dup K (~u) , (39) 



such that r/ + (0) = and with the asymptotic behavior obtained from Eq. (|3"0|) 

V+(y) ~ 7r- lo S^ > y^oo. (40) 

Z7T 

We conclude this subsection by noting that, for d — 2, the statistics of real roots of K n (x) is identical in the 4 
sub-intervals [— oo, —1], [—1, 0], [0, 1] and [1, +oo]. Instead, for d ^ 2, the statistical behavior of real roots of K n (x) 
depend on d in the two inner intervals, while it is identical to the case d = 2 in the two outer ones. In addition, we 
will see below that the polynomials K n (x) © take independent values in these 4 subintervals. 



2. Weyl polynomials 

For Weyl polynomials W n (x) in Eq. ([6]), the expression of the correlation function (|16|) together with the expression 
for the density in Eq. (f2"Tj) yields 



1 / x 2n (x 2 -n- 1) x in + 2 

where T(n,x) = dte~ t t n ~ 1 is the incomplete gamma function [27|]. In Fig. [21 we show a plot of p n {x) (14"T]) for 
different values of n = 50, 100 and 500. One obtains straightforwardly, in the limit n — > oo the uniform density 

Poo Or) = - • (42) 

For n large but finite, the density is uniform like in Eq. (|4"2l up to \x\ ~ y^n above which it vanishes (see Fig. [2]). 
Indeed, one shows in Appendix [Bl that for n ^> 1, one has 

/ s J 71 "" 1 > |ar| 



One notices that this behavior of the density of real roots for Weyl polynomials (|43|) is similar to the density of real 
eigenvalues for Ginibre random matrices (30|, i.e. random nx n matices formed from i.i.d. Gaussian entries. Besides, 
from this scaling form (|43j) one obtains the number of real roots in the interval [—x, x], x > 0, in the large n limit as 

<*„[-*,*]) = f ^ {f j£ ' * <^ (44) 

from which one gets the total number of real roots for n ^ 1 (see also Ref. Slj ) 

2 ^ 

(i\U-oo,+ool) sfn. (45) 

To our knowledge, the higher order terms in this large n expansion are not known. 
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FIG. 2: Mean density of real roots p n {x) given in Eq. (|41|l for Weyl polynomials W n (x) ([6]) as a function of x for different 
values of n — 50, 100 and 500. 

3. Binomial polynomials 

For binomial polynomials, the computation of the density p n {x) is straightforward. Indeed, using Eq. (117|) together 
with the formula for the density (|2ip. one obtains 



p n {x) = Vnp B (x) , p B (x) = 1 

7T(1 + .X z j 



(46) 



exactly for all n > 1 [lj, l2f| . 




FIG. 3: Scaled mean density of real roots n 1 ^ 2 p n (x) as a function of x (|46[) for binomial polynomials B n (x) lff[l. 

In Fig. [31 we show a plot of p B {x) as a function of ir. This formula PB| yields 

(NJa, b}) = [ pjt) dt = ^(ArcTanb - ArcTana) , (47) 

from which one gets the very simple result (see for instance Ref. [13]) 

(A„[-oo,+oo]) = \fn , exactly \fn . (48) 
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B. "Gap probability" on the real axis 

We now study another aspect of the statistical properties of the real roots of these polynomials and focus on the 
probability Po([a, b], n) that they have no real root on a given interval [a, b]. The interval under study will depend on 
the polynomials K n {x), W n (x) or B n {x). 

1. Results from the correlation function 

These polynomials, as a function of x, are Gaussian processes and therefore their zero-crossing properties are 
completely determined by the two-point correlators given in Eq. (|15H17jl . 

Generalized Kac polynomials: For these polynomials K n (x), given the singularity of the mean density p n (x) 
around x — ±1 (see Fig. [T]), it is natural to study separately Po([0, x], n), for x < 1, and Pq([x, oo], n) for x > 1. 
We first focus on Po([0, x], n) and reparametrize K n (x) with a change of variable, x = 1 — 1/t. One finds that the 
relevant scaling limit of C n (t, t') is obtained for t, t' , n — ► oo keeping t = t/n and v = t' jn fixed. In that scaling limit 
the discrete sum in Eq. (115|) can be viewed as a Riemann sum and one finds 

C n (x, y) « n^l m _ x + , (49) 

where I m (y) is defined in Eq. (|29|) . Thus the normalized correlator C n (t, t') — > C(t, t') with the asymptotic behaviors 
(see Eq. §M§ ) 

( d 

Cftivlfe) 4 ' (50) 

t, t' > 1 . 

Thus this correlator is exactly the same as the one found for diffusion, C(t, r) = a(t, t') in Eq. (TTZ)) . Since a Gaussian 
process is completely characterized by its two-point correlator, we conclude that the diffusion process and the random 
polynomial are essentially the same Gaussian process and hence have the same zero crossing properties. Therefore, 
in complete analogy with Eq. ([2]) we propose the scaling form for generalized Kac polynomials 

P (lO,x],n) = A^ n n- m h-(n(l-x)) , (51) 

where A^ n , which is independent of x, is such that linin^oo log„4^" n / logn = and h~{y) — > 1 for i/ « 1 whereas 

h~(y) ~ for !/> 1, where 9(d) is the persistence exponent associated to the diffusion equation in dimension d. 
Defined in this way (fSTj) . h~(u) is a universal function (see below), although the amplitude A~ Ln is not. Note that n 
here plays the role of L 2 in diffusion problem while the variable 1 — x is the analogue of the inverse time 1/t. 

Similarly, we focus on Pq ( [x, oo] ,n), x > 1 , and reparametrize the polynomial with a change of variable, x = 1 + 1/t. 
One finds that the relevant scaling limit of C n (t,t') is obtained for t,t',n — > oo keeping t — t/n and t' = t'/n fixed 
and one obtains 



(52) 



where I m (y) is defined in Eq. (f29|) . Thus C n (x,y) — > C(t,t') with the asymptotic behaviors (see Eq. (|A5j) ) 



C(i,?Hl 4 (TwJ ' M,<<1 (53) 
[l, 

independently of d. Therefore, in complete analogy with Eq. ([2]) we propose the scaling form for random polynomials 

P ([x, oo],n) = ^n-^^+Wz - 1)) , (54) 
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where , which is independent of x is such that limj^oo \ogA\ n / logn = and h + (y) — ► 1 for y -C 1 whereas 
~ y e(2) for y > 1. 

Using the correlator C n (x,y), one shows the statistical independence of the real roots of K n (x) in the four sub- 
intervals [—oo,—l], [—1,0], [0,1] and [l,+oo]. Consider for instance the intervals [0,1] and [l,+oo]. Given that the 
real roots in the interval [0, +oo] are concentrated, for «» 1, around i = lwe introduce x = 1 — 1/ 1 and y = 1 + 1/t' 
and consider the limit t, t' , n — > oo. One easily obtains 

C n (l - 1/t, 1 + 1/t') oc e -«/M«(t,t') ( (55) 

which decays to exponentially for large n. Therefore one concludes that the zeros of K n (x) in the sub-intervals 
[0, 1] and [l,oo] are essentially independent. In a similar way, one shows that the real roots of K n (x) on the four 
subintervals delimited by ±1 are statistically independent. Finally combining Eqs ([5T1 [52} together with ([55| one 
obtains the exact asymptotic result for the probability of no real root as 

F ([-oo,oo],n)ocr^ 2 W d )+ 9 ( 2 » . (56) 

We conclude this paragraph by presenting a heuristic argument which allows to connect the zero crossing properties 
of the diffusion equation to the one of the real roots of K n (x). For that purpose, we consider the solution of the 
diffusion equation with random initial condition (|11[) and we focus on (/>(0, t), without any loss of generality. Following 
Ref. [24|, one observes that the solid angle integration in that expression can be absorbed into a redefinition of the 
random field, yielding 

s l/2 rL 



<K°> *) = I dr r^/ 2 e- r /*¥(r) , (57) 



where Sd is the surface of the <i-dimensional unit sphere and ^S(r) is given by [24 



,. 1 



Ar^O A 



9{r) = S7 1/2 r~^^ Urn -f- / dx ^(x) , (58) 



r<x<r-\-Ar 



which is thus a random Gaussian variable of zero mean and correlations (W(r)W(r')) = S(r — r'). Performing the 
change of variable u = r 2 in Eq. (|57|) . one obtains 

L 2 

<p{0,t)(x[ duvTf L e- u/t ^{u) , (f(it)f(u')) =6(u-u') . (59) 



On the other hand, if we focus on the real zeros of K n (x) in the interval [0, x) with x < 1, we know that these zeros 
accumulate in the vicinity of x — 1. Therefore, in terms ofx = 1 — 1/t one has 

n 

K n (x) ~ a + * (d - 2)/4 ^ l/ *a l . (60) 
»=i 

By approximating the discrete sum in the above expression (|60p by an integral, one sees that K n (x) is similar to the 
solution of the diffusion equation in Eq. (|59[) where L 2 is replaced by n and 1 — x by 1/t. Therefore one understands 
qualitatively why the zero crossing properties of these two processes coincide. 

Weyl polynomials: To analyse the correlation function in Eq. P^)l in the large n limit we write it as 

C n (x,y) = (W n (x)W n (y)) =E^r= ^vf^M^ ' (61) 

' i\ 1 (n + 1) 

i— ' 

where the last equality can easily be obtained using the recursion relation T(n+1, z) = nT(n, z) + e~ z z n . The behavior 
of r(n, z) for large n is analysed in detail in Appendix iBl From the results obtained in Eq. (|B3I IB4[) . one sees that 
the correlation function C n (x,y) in Eq. (|61[) behaves differently for xy < n and xy > n. 

For xy < n, Eq. (|B3[) shows that T(n + 1, xy) — * r(n + 1) for large n so that one finds that C n (x, y) — > C(x, y) with 

C(x,y) = exp[-i(x-y) 2 ] . (62) 
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Interestingly Eq. (f6"2")l shows that inside the interval [—y/n, y/n\, W n (x) is exactly the GSP characterizing the zero 
crossing properties of the diffusion field in the limit of infinite dimension d —> oo (|14p . Therefore one expects 
Po{[— x, x},n), the probability that W n {x) has no real root in the interval [—x, x], with 1 <C x < \fn, to behave as 

Pq{[— x, x], n) (x exp(— 2# oc x) . (63) 

For xy > n, the behavior of C n (x,y) is quite different. Indeed, using the asymptotic behavior in Eq. {B4j, one 
shows that the relevant scaling limit is obtained for x,y,n — > oo keeping x = x/y/n and y — y/\fn fixed such that 
C n (x,y) — > C(x, y) with 

C(i,y) = ^34^ ,£y>l,. (64) 
aty- 1 

Performing the change of variable x — > x = ^/n + X/{\/nt] one easily obtains that C(t,t') behaves like in Eq. {53}. 
Therefore, by analogy with Eq. (|54jl . one deduces that, for < x — y/n « 1, n » 1 keeping y/n(x — \fn) fixed, one 
has 

Pq{\x, oo],n) oc n 2 w(\/n(x — vn)) > (65) 

where w(u) ~ c st for u <C 1 and tu(u) cx u 6 ^ 2 ' for «> 1. In addition, following the arguments presented above (see 
Eq. I|55p). one shows that these two outer-intervals [—00, —\/n\ and [y/n, +00] are statistically independent, such that 

P ([-oo, — v/n] U [Vn,+oo],n) oc n^ 0{2) . (66) 

However, given the behavior of the correlator C(x,y) for xy > n in Eq. (|62p the inner and outer intervals are not 
independent and the probability of no root on the real axis is not the product of the probabilities in Eq. (|63p evaluated 
in x — \fn and the the one in Eq. ((66)) : the effect of these correlations will be discussed below. 

Binomial polynomials: In that case, one can extract information directly from the correlation function in Eq. (fT7|) 
by focusing in the limit x, y — > 0. In that limit the normalized correlation function C n (x, y) is given by 

C n (x,y) = + n ~exp[-§(x-y) 2 ] ,x,y<Zl. (67) 

[(l + z 2 )(l +2 / 2 )]2 

Thus in the large n limit, the probability Po([a, b], n) that binomial polynomials have no real root in the interval [a, b] 
1 

with a<k<l, n 2 «i-s behaves like 

Po([a, b],ri) oc exp [-0ooVn(b -a)] , a < b < 1 , n~2 < & - a , (68) 

which is an exact statement. However it is a more difficult task to obtain the behavior of Po([a, b], n) for an arbitrary 
interval [a, b] and eventually obtain the probability of no root on the entire real axis for this class of polynomials ([7]) : 
this will be achieved in the next sections. 

To conclude this paragraph, we have shown that the analysis of the correlation function C n {x, y) yields important 
exact results for the gap probabilities. Indeed, for generalized Kac polynomials, we obtained the important results in 
Eq. ([51]) and Eq. (|54|) which yield the exact result in Eq. (|56|) . For Weyl polynomials the study of the correlation 
function allowed us to obtain the results in Eq. (|63j) and Eq. (|65[) . Finally, for Binomial polynomials we obtained, 
from the correlation funtion, the asymptotic behavior in Eq. (|68p . which will be useful in the following. 

2. Mean-Field description : Poisson approximation 

To calculate the gap probabilities and the associated scaling functions, we first develop a very simple mean field 
theory. This theory, albeit approximate as it neglects the correlations between zeros, is simple, intuitive and qualita- 
tively correct. We will see later how one can improve systematically this mean field theory to get answers that are 
even quantitatively accurate. As a first step, we neglect the correlations between the real roots and simply consider 
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that these roots are randomly and independently distributed on the real axis with some local density p n (x) at point 
x. Within this approximation the probability Pk{[a,b],n) that these polynomials have exactly k real roots satisfies 
the equation 

dP k+1 ([a, b],n) = pn{b)[ p k{[ai b]f n) _ p k+i{[a> b]>n)] > (69) 

together with the normalization condition X)fc>o Pk([a, b], n) = 1 and Pk([a,b],n) — 5k,o when a = b. In the large 
n limit (where one can omit the constraint Pi- >n ([a, b], n) = 0), Pk{[a,b],n) is given by a non-homogeneous Poisson 
distribution 

P k ([a,b],n) = j^e-" , n={N n [a,b]) = J p n {x)dx , (70) 
which clearly satisfies Eq. (l69|) . In particular, this mean field approximation (|70|) yields the gap probability 

P ([a,b],n) = exp(-(N n [a,b])) = exp \- J p n (x) dx^j . (71) 

When applied to Generalized Kac polynomials K n (x), for which we obtained (N n [0, x]) in Eq. (|35[) . this mean-field 
approximation (|7ip yields in the scaling limit n — > oo, 1 — x — > 0, n(l — x) > fixed 

_j_ fd 

P o ([0,x],n)=Al n n 2 - V 2 exp [77.(71(1 - x))} , (72) 

where, from Eq. (f3"2"|) . log^l^^ = o(logn). This mean-field approximation thus yields the correct scaling from for 
Po([0, x], n) as in Eq. (I51j) . with the non trivial predictions for the exponent and scaling function 



d2/P A (y)J, (73) 
with the asymptotic behavior h~(u) ~ 1 for u <C 1 and, using the asymptotic behavior obtained in Eq. (|37ll . 

M M )~ u * MF ( d ). 

Similarly this mean-field approximation applied to Weyl polynomials W n (x), for which we obtained (N n [—x, x]) in 
Eq. (jH}, yields the scaling form 

M-*.*],n)~Hf-f^'*<^L (74) 
I exp (—Zy'n/ir) , a; > ^Jn . 

Notice that this mean-field approximation gives an approximation of 6*^ F = ir^ 1 = 0.31831... (see Eq. (|6"3")0 . which is 
consistent, using the relation 9(d) = 2~ z / 2 6 O0 ^fd for large d, with Eq. ([73]) . 

Finally, if one uses this mean-field approximation to study binomial polynomials B n (x), one obains, using the 
expression (N n [a, b]) given in Eq. I|47p 



P ([a,b],n) = exp 



- — ( ArcTan a - ArcTan b) 

7T 



which again, according to Eq. ([68]) . gives the mean-field approximation for 8^ F = n 1 , as above. 



(75) 



3. Beyond Mean-Field : a systematic approach 

We will now show that this mean field approximation (|7ip can actually be improved systematically. For that 
purpose, one considers the probability Pk(\a, b],n) that such polynomials as in Eq. ([3l[6j[7]) have exactly k real roots 
in the interval [a, 6]. Following Ref. [21(, one introduces the generating function 

00 

P n (p,[a,b]) = J2P kp ^ a ' b ^ n ^ ( 76 ) 

fc=0 
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where P n {p, [a, b]) can be interpreted as a persistence probability with partial survival [2l| . For a smooth process, it 
turns out that 9 n (p, [a, b]) = — log (P n (j>, [a, b])) depends continuously on p : this was shown exactly for the random 
acceleration process (see Eq. (|106p below) and approximately using the IIA - and further checked numerically - for 
the diffusion equation with random initial conditions [2lj . Thus one has 

§ n (p, [a, b}) = log (P n (p, [a, b})) = £ (l0g ( f )r (iV;([a, b])) c , (77) 

r=l 

where the notation (...) c stands for a connected average. Here we are interested in P n (p = 0, [a, b]) = Po([a, b] 7 n) and 
the idea, given that 9 n (p = 1, [a, b]) = is to expand n (p, [a, b\) around p = 1 in an e-expansion with p = 1 — e. This 
yields 

oo 

§ n (l-e,[a,b])=J2*r,n(fab]y , (78) 

where a rjTl ([a, 6]) are linear combinations of the cumulants (iV™[a,6]) c , with m <r. For instance, 

CH,„(MD = (N n [a,b\) , a 2 ,„([M]) = \ {{N n [a,b]) - Ma,&]) c ) , (79) 
a3,n([o, 6]) - 3 g + g 

Thus one sees that if one restricts the e-expansion in Eq. (|78|) to first order, and set e = 1, one recovers the mean- 
field approximation (I7ip . using Oi jM ([a, 6]) = (iV n [a, 6]). Higher order terms in this e-expansion allow to improve 
systematically this mean-field approach. 

Kac polynomials for d = 2: We first illustrate this e expansion for Kac polynomials K n {x) for d = 2 where we 
compute 9 n (l — e, [0,x]) up to order 0(e 2 ). In that purpose, we compute a2, n ([0,a;]). In Appendix [Cl we show that 
in the scaling limit (1 — x) — > 0, and n — > oo keeping the product n(l — x) fixed one has, similarly to the Eq. (|35[) for 
the first moment, 

(JV2[0,x]) c = <A^[0, l]) c - u_(n(l - x)) (80) 
where V- (y) , given in Eq. (|C14|1 , is such that (y) for y <C 1 and 

v-{v)~ (~--y\logy, y»l. (81) 

Notice that (N%[0, l]) c in Eq. ([80]) has been computed in Ref. [32|, yielding for large n 

(^[0,l]>c= (--^)\ogn + o(logn) , (82) 

although higher order terms in this large n expansion are not known. Combining Eq. (|78l [79l [80]) together with the 
expression for v~(y) in Eq. (IC14j) . one obtains in the scaling limit 

§ n (l-e,[0,x])= (e+4) <A^„([0,1]))-4(7V„([0,1]) 2 ) C (83) 
-e [ dup K (u) + ( 2 f dux ( du 2 (}C(u 1 ,u 2 )- p K (ui)p K (u 2 )) +0(e 3 ) , 

JO JO Ju t ^ ' 

with y — n(l — x). In this Eq. (|83[) . p K (u) is given in Eq. (|3ip . and /C(ui, 112), which we study in detail in Appendix [Cl 
is essentially the two-point correlation function of real roots inside the peak of the density around x — 1 (see right 
panel of Fig. [T|). Setting e = 1 in this expression to order 0(e 2 ) (|53")) one obtains, for d = 2 

P ([0,x],n) =A- d n n- e{2) h-{n{x-l)) , 0{2) = = 0.180899... , (84) 

h-{y) = exp(^J\up K {u)- J" d Ul J du 2 {t{ui,u 2 )- p K {u 1 )p K {u 2 )) S j , (85) 
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with \ogA d n = o(logn) (see Eqs ([521 |55J|). Note that the value of the exponent 8(2) up to second order as given in 
this Eq. (I84p was computed in Ref. [21]. In the next sections, we will show, using numerics, that this second order 
calculation (|84p is a true improvement upon the mean-field approximation (|73[) . 

Weyl polynomials: We now compute 8 n (l — e, [—or, x]), x < y/n for Weyl polynomials up to order 0(e 2 ). In that 
purpose we compute fl2,n([ — x, x\). In Appendix iDl one shows that for x < y/n fixed and n ^> 1, one has 

(N*[-x,x}) c = v(x) , (86) 

with v(x) given in Eq. (|D1 1|1 . Combining Eq. ([751 1751 [56]) together with the expression for v(x) in Eq. (|Dlip . one 

obtains for large n 

6Jl - e, \-x, x]) = e— + 2e 2 / ds(s - x)(W(s) - ^ 2 ) + 0(e 3 ) , (87) 

I" JO 

where W(s), given in Eq. (|D8p . is the two point correlation function of real roots in the interval [— y/n, y/n}. Setting 
e = 1 in this expression up to order 0(e 2 ) ([87]) . one obtains, for x < y/n 

P ([-x,x],n) =exp (~ -J ds(s - 2x)(W(s) - tt" 2 )^ . (88) 

Using Eq. (|D13p . one obtains its large x behavior as in Eq. ([6"3"]) with the value of 8^ up to order 0(e 2 ) 

1 C°° 

6oo = -- \ ds(W(s) - it- 2 ) = 0.386471... , (89) 
Jo 

which should be compared with the numerical value ^oo,sim = 0.417(3) 

If we focus instead on the outer intervals [— oo, —y/n\ U [y/n, oo], this e expansion is essentially similar to the one 
performed for Kac polynomials and d = 2, given the behavior of the correlator in Eq. (|64|) . More interestingly if we 
are interested in the computation of the probability of no root on the real axis, this e expansion is able to take into 
account (perturbatively) the correlations between the inner and outer intervals, which, as discussed below Eq. (|66l) . 
can be seen in the correlation function ([6"2"| . Doing so, one obtains that 

P ([-oo,oo],n) = P Q ([-y/n, y/n],n)P ([-co, -y/n] U [y/n,oo],n)p n , (90) 



where the computation of p n is similar to the one carried out in Appendix [Dl where the quantity W„(ii, t^) in Eq. (|E)2|) 
involves t\ € [—y/n, y/n\ and £2 € [—00, — y/n] U [y/n, 00} such that £1^2 < 1. The computation up to order 0(e 2 ) shows 
that p n oc n T with r > 0. Therefore, on the basis of this result together with Eqs ([63] 166 ]) . one expects 

P ([-oo, 00], n) ~ n~ 7 exp (-26^^) , (91) 

where the exponent 7 is a priori unknown. Below, we will confront this statement with numerical simulations. 

Binomial polynomials: We now focus on 8 n (l — e,n) for binomial polynomials. In Ref. [20j, it was shown that 

1 

for large n and all a,b > n 2 

(iV 2 [a,6]) c cx/3 2 (7V„[a,6]) , (92) 

where (3% is a constant, independent of n, a and b. It has an analytic expression in term of an integral 
involving elementary functions, with /3 2 = 0.571731.... This expression (|92[) yields in the large n limit, 
&2,ra([a, b}) = (1/2 — /3 2 ) (N n [a, b]) which together with the asymptotic behavior in Eq. ([68]) and the cumu- 
lant expansion of Eq. ([78]) allows to compute 9^ up to order 0(e 2 ). We have checked that this coincides with the 
one obtained in Eq. (155]) . More generally, one expects that for all integer m > 

(N™[a,b}) c <xp m (N n [a,b}) , (93) 

where (3 m is a constant, independent of n, a and b and in Appendix [El we explicitly show the mechanism leading to 
this relation (|93p for k — 3. From these relations (]Mj) and the structure of the cumulant expansion in Eq. ([75p . one 
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expects that Po([a, b],n) oc exp [— u>(N n [a, b])], where a; is a linear combination of the coefficients (3 m . Finally, this 

_ l 

expression has to match the exact asymptotic behavior of Po([a,b],ri) for a < b <C 1 and n 2 <C b — a derived in 
Eq. (j6"5)) . Thus one has u> = irOoo so that one obtains the exact result 

Po([a, b],n) oc exp [-^/n9 00 (ArcTanb - ArcTana)] , (94) 

from which we obtain the exact expression for the probability of no real root for B n (x) in the large n limit as 

Pq([— oo, +oo], n) oc exp (— v^tt^oo) • (95) 

Below, we check this analytical result (194)) using numerical simulations. 

4- Numerical results 

Kac polynomials: We first focus on the interval [0,1] and check numerically the scaling forms for Po([0, x],n) 
in Eq. ([2]) for K n (x) and for different values of d. In each case, this probability is obtained by averaging over 
10 4 realizations of the random variables a,i 's in Eq. ([3|) , drawn independently from a Gaussian distribution of unit 
variance. In Ref. Jjj], we already presented numerical results for Pq([0, x],n) in d = 2 and we also checked that 
P ([0,x],ri) oc n~ 6 ^ with 0(2) = 0.187(1). In the left panel of Fig. gj we show a plot of log[P ([0, x], n)/P o ([0, 1], n)] for 
d = 2 as a function of the scaled variable n(l — x) for different values of n. According to Eq. |2]), together with the good 
collapse of the curves for different values of n, this allows for a numerical computation of the scaling function h~ (y). 
We have checked that different distributions of the random coefficients either a% — ±1 or rectangular distribution yield 
the same scaling function h~ (y), suggesting that this function is indeed universal. On the same figure, the dotted line is 
the result of Mean Field approximation ([751) i or nrs t order in the e expansion, and the solid line is the analytical result 
of the second order calculation obtained in Eq. (|84|) . In both cases, the integrals involved were evaluated numerically 
using the Mathematica. As expected one observes on this plot that the Mean Field calculation is only in qualitative 
agreement with the numerical results, we recall in particular that #mf(2) = 1/2-7T = 0.159155.... Interestingly, one 
sees that the second order calculation is a clear improvement over the Mean Field calculation which is in quite good 
agreement with the numerical results for the scaling function, in particular, 9(2) — (tt + 4) / Air 2 = 0.180899.... We have 
checked that this scaling ^ holds for other values of d. In the right panel of Fig. [U we show a plot of n e ^ 3 'Po([0, x],n) 
for d — 3 and 9(3) = 0.238(4) as a function of n(l — x) for different degrees n — 256, 512, 1024. Again, the value of 
9(3) = 0.238(4) for which one obtains the best collapse of the curves for different values of n is in good agreement 
with the values of 9(3) found for the diffusion equation [1, 0|- 

We have also checked numerically our results for the gap probability in the outer intervals (|54p . For that purpose, 
we notice that Pq([x, oo], n) — Po([Q, 1/x], n), which is easier to compute numerically, where Po([a, b], n) is the gap 
probability associated to the polynomial K n (x) defined such that K n (x) = x n K n (l/x) with 

K„(x) = ^2a n -i(n-i) 4 x l +a x n . (96) 

i=0 

Thus, from Eq. ([54]) one expects that for < 1 — x <C 1, and n ^> 1, keeping the product n(l — x) fixed one has 

P Q ([0,x],n) (xn- 0{2) h + (n(l-x)) , (97) 

independently of d. In Fig. O we show a plot of n 9 ^ P a ([0, x], n) as a function of n(l — x) for d = 3 and for different 
values of n = 128, 256, 512, 1024. Again, the good collapse obtained for 9(2) = 0.1875 corroborates the validity of the 
scaling in Eqs (HI E7J- 

Weyl polynomials: We first focus on the inner interval [—^Jn, y/n] and compute numerically the gap probability 
Po([— x, x], n) for x < yjn. Here also, this probability was computed by averaging over 10 4 different realizations of 
the random coefficients <2j's. In the left panel of Fig. [HI we show a plot of [log Po([—x, x], n)]n~~ 1 / 2 as a function of 
xj\Jn < 1 for different values of n = 40,90 and 150. According to our prediction in Eq. ([63| . Po([—x,x],n) behaves 
exponentially for large x. From the slope of the straight line in the left panel of Fig. one extracts 29^ = 0.845(3), 
in good agreement with previous numerical estimates from the persistence probability for the diffusion equation in 
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FIG. 4: Left : Plot of log[Po([0, x],n) /Po([0, 1], n)] for Kac polynomials K n (x) and d = 2 as a function of n(l — x) for different 
degrees n = 256, 512, 1024. The dotted is the result of the Mean Field prediction (|73[) and the solid line the result of the second 
order calculation in the e expansion ([83]). Right : Plot of n e(3) P o ([0, x],n) on a log- log scale for Kac polynomials K n (x) and 
d = 3 with 6>(3) = 0.238(4) as a function of n(l - x) for different degrees n = 256, 512, 1024. 
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FIG. 5: Plot of n e< - 2 'Po([0, x], n) on a log-log scale, where K n (x) is defined in Eq. (|96|l as a function of n(l — i) for d = 3 and 
for different values of n = 128, 256, 512, 1024. 



large dimension [6j, |7j. On the same figure, we have also plotted with a dotted line the result from the Mean Field 
approximation (|74|) and with a solid line the result up to second order in the e expansion in Eq. (|87|), Again, the 
second order term allows to improve significantly the Mean Field prediction. We recall the estimate up to order C(e 2 ), 
26*oo = 0.7729... from Eq. (f89|) . We now focus on the outer intervals [— oo, — y/n] and [y/n, +oo] . In the right panel 
of Fig. we plot n e ( 2 )Po([— x, — \fn] U [\fn, x},n) as a function of x/^/n > 1 for different degrees n = 20,40 and 
90. The fact that the curves for different n collapse on a single master curve is in agreement with the scaling 
proposed in Eq. (prj)) . Finally, we computed the gap probability on the full real axis. In the right panel Fig. [71 
we plot [log(n 7 Po([— x, x], n))]n -1 / 2 as a function of x/^/n for different values of n — 40,90 and 150. The exponent 
7 = 0.10(1), which is the only fitting parameter is fixed to obtained the best collapse of the different curves in the large 
xj\fn limit. The solid line has a slope —28^ = 0.845, which is also the value reached by log(n 7 Po([ — x, x], n))n~ 1 / 2 
for large x. This fact is in complete agreement with the scaling proposed in Eq. (f9"Tj) . The fact that 7 < 9(2) arises 
from the correlations between the inner and outer intervals. 

Binomial polynomials: Finally, we have checked the exact result in Eq. (|94|) for binomial polynomials ([7]). For 
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FIG. 6: Left : Plot of [log Po([ — x , x], w)]n _1,/2 for Weyl polynomials W„(x) ([6]) as a function of x/^/n < 1 for different values 
of n — 40, 90 and 150. The dotted line is the result of the Mean Field approximation whereas the solid one is the result of 
the expansion up to order 0(e 2 ) in Eq. (57) . Right : Plot of n e ^ 2 'Po([— x, — y^n] U [\/n, x],n) as a function of x/y/n > 1 for 
different degrees n = 20, 40 and 90. 
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FIG. 7: Left : Plot of [log(n 7 P ([-z, x], n))]n~ 1/2 for Weyl polynomials W n (x) © as a function of x/ y/n for different values 
of n = 40,90 and 150. Here 7 = 0.10(1) and the solid line has slope -0.845. Right : Plot of [log P ([-x, x], n)]n~ 1/2 for 
binomial polynomials B n (x) ((Jj as a function of x for different values of n = 40 and 90. The solid line is the exact result in 
Eq. (I94|l with #oo = 0.42(1), which is the only fitting parameter here. 



that purpose, we have computed numerically Po([— x, x], n) by averaging over 10 different realizations of the random 
coefficients dj's. In the right panel of Fig. [Jj we show a plot of [logPo([— x, x], n)]n~ 1 / 2 as a function of x for different 
values of n = 40 and 90. The solid line is the analytic prediction from Eq. with = 0.42(1), consistent with our 
previous estimates from Weyl polynomials (see left panel of Fig. [7j. The good agreement with the numerics confirms 
also the exact result for the probability that these polynomials have no real root on the real axis in Eq. (|55|) . 



C. Probability of k real roots : large deviation function 



We now generalize our analysis and study the probability Pk([a,b],n) that such polynomials El EJ) have exactly 
k real roots [i~8| in a given interval [a, b]. 
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1. Mean field approximation : Poisson approximation 

One first considers the Mean Field approximation introduced above where one assumes that the real roots are 
totally independent and randomly distributed with density p n {x). This leads to Eq. (|70[) 

ft([a,6],n)= <jVm[ °,' 6])fc e-<W>. (98) 
fc! 

If we focus on the limit n ^S> k ^S> 1, keeping the ratio y = kj (N n [a, b]) fixed, one has 

log P k {[a, b],n) (N n [a,b]) ip M F [ ( N [ a ty j > VMf(j/) = 1 + V logy - y , (99) 

where we have used the Stirling's formula log (fc!) = (k+ 1) logfc — k + 0(l). We will see below (through scaling analysis 
as well as numerics) that this Mean Field approximation provides the correct scaling form for P k ([a,b],n) (although 
the exact computation of (p(y) certainly demands a more sophisticated analysis). Let us present the consequences of 
this scaling form in Eq. (|99p for the different polynomials under study. 

Kac polynomials: Let us define qk(n) — -Pfe([0, 1],ti). In that case, we have seen in Eq. (|32|) that 

(iV n ([0, 1])) ~ ^:^/| log n so that one expects the scaling form 

( [2 k \ 

log q k {n) (x -log nip \2tv\ -- . (100) 

\ \ a log n I 

For the special case of Kac polynomials (d = 2), this scaling form, in the neighborhood of fc = logn/27r, is consis- 
tent with the rigorous result [321 ] that in this neighborhood q k (n) is a Gaussian with mean logrt/27r and variance 
V n ~ - Jr)logn ([82]) in the large n limit. 

Weyl's polynomials: Let us define qk(ri) — Pfe([— oo, oo], n). According to Eq. which tells us that 

(N n [— oo, oo]) ~ (2/7r)i/n, and the scaling form in Eq. ([99]) one expects 

logg fc (n) oc -«/n<p [^^=j ■ ( 101 ) 

Binomial polynomials: Let us define similarly qk{n) — Pfe([— oo, +oo], n). In that case, according to Eq. ([48]) 
which tells us that (N n [— oo, oo]) = y/n, and the scaling form in Eq. (l99|) one expects 

logq k (n) oc -y/nip (~^J ■ ( 102 ) 
In the following, we will check these scaling forms P00H102")) numerically. 



2. A more rigorous approach for a smooth Gaussian stationary process 



We illustrate this approach on the diffusion equation with random initial conditions ([T]), which is the underlying 
stochastic process describing the statistics of real roots of these random polynomials. We thus consider the probability 
Pk(t,L) that the diffusing field </>(x, t) crosses zero exactly k times up to time t. Let us first consider the regime 
1 t <C L 2 . In this regime, p k (t,L) is given by the probability Vk(T) that X(T) crosses zero exactly k times 
where X{T) is a GSP with correlations a(\T - T'\) = [cosh(|T - T'|/2)]- d / 2 , where T = logt. Since, a(T) = 
1 — jqT 2 + o(T 2 ) for small T, this GSP is a smooth process with a finite density of zero crossings given by the Rice's 
formula /i — [— a" (O)] 1 ^ 2 /n 26]. We propose the following scaling form for large T and large fc, with k/T fixed 



logP fc (T) = - 




(103) 
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To understand the origin of this scaling form, let us consider the generating function V{p, T) = X)a*1o P^C ' T) as in 
Eq. ((76|) . One can show 21] that V(p, T) ~ exp[— 8(p)T], where for a smooth GSP 0(p) depends continuously on p. If 
the scaling in Eq. (|103[) holds, one gets by steepest descent method valid for large T, 8(p) = Min x>0 [[ix \ogp — <p(x)]. 
Inverting the Legendre transform we get 

<p(x) = Max < p < 2 [^logp + §(p)\ . (104) 

Notice that although 6(p) is a priori defined on the interval [0,1], the computation of (p(x) involves an analytical 
continuation of 8(p) on [0, 2). Going back to real time f, Eq. (| 1 03[) then yields a rather unusual scaling form valid in 
the limit 1 « t « L 2 

logp^L^-logf^^) . (105) 

In the opposite limit t^>L 2 , one simply replaces f in (|105|) by L 2 . Translating into random polynomials, this regime 
corresponds to (1 — x) <C n~ 1 since one just replaces f by 1/(1 — x) and L 2 by the degree n as discussed before. 
Hence, in this regime, we arrive at the announced scaling form for qk{n) in Eq. (jlOOp . This approach can be extended 
straightforwardly to the other classes of polynomials, yielding the scaling forms in Eq. (|10H I102|) . 

Of course, despite the exact formula (| 1 04[) . the function <p(x) remains very hard to compute, simply because 9(p) 
is, in many cases, unknown. However, for the random acceleration process (RAP), sometimes called in the literature 
"integrated brownian motion" 

£jjW=r,(t), ( V (t) V (t'))=S(t-t'), (106) 

where 77(f) is a white noise for which /xrap = %/3/ (27r), #rap(p) has been computed exactly [H,[34|], yielding #rap(p) = 
i[l — ~ sin (§)]. By performing the Legendre transform (|104|) one obtains 

wapW - fofes) + 1 - ; AreSi " taffs)) • <107) 



with the asymptotic behaviors 



^rap(*)~ i ll^-l) 2 , x^l (108) 



4 , 2 J-xlogx , x —> 



^-x log 2 , x — > 00 
which gives back the exact result </?rap(0) = 1/4 (35l. [36|. 



3. Numerical results 

In Ref. [TH, we have checked numerically the scaling form (|105[) for the diffusion equation with random initial 
conditions. Here, we have computed numerically these probabilities qk(n) for the different polynomials under study. 
This was done by averaging over 10 4 different realizations of the random coefficients a^s. In the left panel of Fig. [8l 
we show a plot of — [log qk(n)]/ \ogn as a function of 27rfc/logn for Kac polynomials ([3]) with d = 2 and for different 
values of n = 20,40,80 and 100. This suggests that the different points fall on a single master curve, which is in 
rather good agreement with the scaling form proposed in Eq. (|100|) . Similarly, in the right panel of Fig. we show 
a plot of — [log qk{n)] rT 1 ! 2 as a function of k/y/n for binomial polynomials ([7|) for different values of n = 31,63, 127 
and 255. Here also, the fact the points fall a single master curve is in good agreement with the scaling form in Eq. 
(|102| . We have also checked that a similar scaling, as in Eq. (jlOlj) holds for Weyl polynomials ([6]). 



D. Distribution of the maximum of real roots 



Up to now, we have mainly focused on the distribution of the minimum A m i n of the absolute values {|Ai|} of the 
real roots of these polynomials. Indeed, the gap probability Po([— x, x], n) is just the probability that A m i n is larger 
than x > 0. We now instead focus on the distribution of the maximum A max of the {|Ai|}. 
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1, Mean Field approximation 

As a first approach, we consider the Mean Field or Poissonian approximation introduced above, where one neglects 
the correlations between the real roots. Then, for any random polynomial Q n (x) = X^o^ 2 ^' ^ ne probability that 
A max < x is simply the probability that Q n (x) has no real root in [— oo, —x] U [x, +oo]. Therefore, within the Mean 
Field approximation one has 



Proba. (A max < x,n) = exp ^-2 J p n (t)dt\ . 



(109) 



Taking the derivative of this expression (| 1 09|) with respect to x, one obtains the probability distribution function (pdf) 
Pma.x(x,n) of the maximum A max as 



Pmax(a;, n) = 2p„(x)exp y~ 2 J Pn{t)dt J . (110) 
To obtain the large x behavior of p m:ix (x,n) one computes the one of p n (x) whose expression is given by 



Pn (x) = ^)(w+y-Kwi a ( Cn{x) = ^ x2i 

ZT7C n (X) 



For large x, one has c n (x) = x 2n (b 2 n ) + x 2n - 2 {b 2 n _ x ) + 0(x 2n ~ 4 ) which gives p n (x) ~ ^ (b 2 ^) / (b 2 ,) / (ttx 2 ). Finally, 
from Eq. (|110p . one obtains for x ^S> 1 



Pu^(x.n) - -IpnU) - ^J ^" ^ . (112) 



irx 



2. Exact result for the tail 

In fact, the tail of the distribution can be computed exactly by noting that 

1 1 



Proba. (A max <x,n) = P 



(113) 



22 



where Po([a, b],n) is the gap probability associated to the polynomial Q n (x) defined such that Q n (x) — x n Q n (l/x). 
Similarly, we denote p n (x) the mean density of real root associated to Q n (x). From this exact identity (|113[) . valid 
for all polynomials Q n {x) and all n > 1, one obtains the asymptotic behavior 

Proba.(A max <x,n) = l- 2^^ + o(x~ 1 ) , (114) 

x 

where we have simply used the definition of p„(0), provided it is well defined, which is the case for Gaussian random 
polynomials. From this asymptotic behavior (|114p , one obtains the pdf p m nx(x,n) of the maximum A max for large 
x > as 



p max (x,n) ~ 2^1 , p n (0) = \\j^y (H5) 

where we have used the formula to compute ~p n (0). For the various polynomials under consideration here, one 

thus obtains such a power law tail (|115p with p n (0) — tt ((n — l)/n) 4 for Kac polynomials, p n (0) = yfnj'n for Weyl 
polynomials and binomial polynomials. This result shows in particular that the mean value of A max is not defined for 
these polynomials. 

It is interesting to note that the Mean Field approximation (| 1 1 2[) gives the exact result for this algebraic tail (| 1 1 5[) . 
This might be understood by noting that, apart from a short range repulsion, the real roots of these polynomials are 
actually weakly correlated. For instance, for Weyl polynomials, we show in Appendix [D] see Eq. (|D10[) , that the two- 
point (connected) correlation function of the real roots inside the interval [—\fri, y/n\ decays faster than exponentially 
at large distance (see also Ref. 20] for similar properties of K n (x) and B n (x)). On the other hand, the marginal 
distribution of a single real root is nothing else but p n (x)/(N n [— oo, oo]) which decays algebraically p n (x) ~ 1/x 2 for 
large x. Given that these real roots are weakly correlated, we expect that the distribution of the maximum of these 
real roots is given by a Frechet distribution, which indeed has a power law tail, as we found here (|115[) . 

3. Numerical results 

We have checked numerically this exact asymptotic behavior (|115[) for the three classes of random polynomials ([3j 
ElE]). In all the three cases the pdf of the maximum A max was obtained by averaging over 10 5 different realizations 
of the random Gaussian variables a^s. In the left panel of Fig. [9l we show a plot of p m ax(^, n) as a function of x for 
Kac polynomials and d = 2, for different values of n = 8, 16 and 32. Notice that the exact result in Eq. (|115[) . which 
is plotted here with a dotted line, is in principle true for all n > 1 so that it is not necessary to perform numerics 
for polynomials of large degree. This figure shows a good agreement between our analytic result and the numerics. 
Similarly, in the right panel of Fig. [HI we have plotted PmaK(x,n)/-\/n f° r Weyl polynomials W n (x) and for different 
values of n = 8, 16 and 32. Here again, the agreement with the analytical result in Eq. (|115p is quite good. 

IV. CONCLUSIONS AND PERSPECTIVES 

In conclusion, we have investigated different statistical properties of the real roots of three different types of real 
random polynomials ((3j[6j[7|). In these three classes, the mean density of real roots exhibit a rich variety of behaviors, 
as shown in Figs [TJ3J We have first focused on gap probabilities (1551 ltJTl I55|) which were shown to be closely related 
to the persistence probability for the diffusion equation with random initial conditions in dimension cZ <jTJ [^]) . We 
proposed a Mean Field approximation to compute the exponents as well as the universal scaling functions describing 
these gap probabilities. Furthermore, we showed how to improve systematically this MF approximation (see Fig. 2J 
|6]) using an e-expansion based on the so called persistence with partial survival. In the case of binomial polynomials 
B n (x) ([7]), this allows to obtain exact results for the gap probability (fM|) . Our main results for the gap probability 
qo(n) on the full real axis are summarized in Fig. 1101 We hope that these connections may help to obtain exact results 
for the exponent 0(d), which remains a challenge. Besides, we extended our analysis to the probability that these 
random polynomials have exactly k real roots in a given interval [a, b\. We have shown, in the three cases, that this 
probability has an interesting scaling form characterized by a large deviation function (|100j|102"]) . We proposed a Mean 
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FIG. 9: Left : Pmax(x,n) for Kac polynomials K n (x) ([3]) and d — 2 as a function of x for different values of n = 8, 16 and 32. 
The dotted line is the exact result for the tail of the distribution (|115[) . Right : p ma x(x, n)/^/n for Weyl polynomials W„(x) 
([6]) as a function of x for different values of n = 8, 16 and 32. The dotted line is the exact result for the tail of the distribution 
(fTT5l) 
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FIG. 10: Summary of the main results for the probability of no root on the full real axis qo(n) for the three different classes of 
polynomials K n (x),W n (x) and B n (x). 



Field approximation which reproduces qualitatively these scaling forms, which were further checked numerically (see 
Fig. [5]). A similar question was asked in the past for Ginibre real matrices : what is the probability that such n x n 
matrices has exactly k real roots ? Quite recently in Ref. [371 ] . Akemann and Kanzieper obtained an exact formula 
for this probability. In that case, the mean number of real roots grows like *Jn and our MF approximation would 
therefore suggest a scaling form for this probability similar to Eq. p01[) . It would be very interesting to extract the 
large deviation function from an asymptotic analysis of the exact result of Ref. [37| . 

Finally, we computed the pdf of the largest real root of these random polynomials. We showed that for a wide class 
of random polynomials, this pdf has an algebraic tail with exponent —2 as shown in Eq. (|115p and it would certainly 
be interesting to extend these results to the case of non-Gaussian random polynomials. 



Acknowledgments 



We thank P. Krapivsky for stimulating discussions. One of us (G.S.) would like to thank the physics department 
of the Saarland University for hospitality and generous allocation of computer time. 



24 



APPENDIX A: CALCULATION OF THE MEAN DENSITY FOR GENERALIZED KAC'S 

POLYNOMIALS 

In this appendix, we give some details concerning the computation of the mean density and the mean number of 
real roots for Kac's polynomials K n {x) 

1. Scaling form 

The starting point of the calculation is the formula for the density p n {x) given in Eq. (|21[) . which using Eq. ([1 
gives 



Pn[ X ) = -J , 8d/ 2 2+UX) f , - ( 1 W 1 S , (AD 

7T Y l + x 2 s d/2 _ l!n (x) \l + a d/2 -i(x)J 



with Sk, n (x) =y^j k 



1=1 



In the limit n» 1, and |1 — x\ -C 1 keeping y = n(l — x) fixed, Sk, n (x) can be viewed as a Riemann sum, thus 



Sk,n{x)~n k+1 [ dxx k e- 2ux , u = n(l - x) . (A2) 
Jo 



Finally using Eq. (|A2|) together with Eq. (|Aip yield the formula (|28[) given in the text : 



Im{y) = I dx x m exp (-2yx) . 
Jo 



(A4) 



2. Asymptotic expansions 



To compute the asymptotic behaviors of the function p K (y) in Eq. (|28[> . one needs the asymptotic behaviors of 
Im(y) defined above (IA4jl . 

a. The limit y — > oo 

In that limit, one writes I m {y) as 

Irn(y) = y' 1 -" 1 duu m e- 2u = y- 1 - m duu m e~ 2u + 0(y m e~ 2v ) 
Jo Jo 

= (2yy 1 ~ m T(l + m)+0(y m e- 2 y) , (A5) 
which, together with Eq. (|28[). gives the first line of the formula given in Eq. (130)) in the text. 
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b. The limit y —* — oo 

Because of several cancellations occurring in p K (y) in the limit y — » — oo one needs the three first terms in the 
asymptotic behavior of I m (y) in that limit : 

f \v\ 

I m (y) = / duu m e 2u 



which after a tedious but straightforward calculations, using Eq. (128|) . gives the second line of the formula given in 
Eq. (1201) in the text. 



APPENDIX B: CALCULATION OF THE MEAN DENSITY FOR REAL WEYL'S POLYNOMIALS 

In this section, we present the details of the calculation leading to the scaling form given in Eq. (|43[) . We start 
with the expression for the mean density given in Eq. (|4ip with x = u^fri: 



. . , _ . 1 / n n+1 U 2n (u 2 - 1 - n- 1 ) n 2n+l u in+2 

Pn{x) = Pn {^u) = -Jl+ enu , T{n+ ^ nu2) - pf(^M^ • (B1) 



We want to obtain the limit n — > oo, keeping u fixed, in the above equation (|B1[) . For that purpose, we need the 
large n behavior of T(n, nu 2 ) for large n and u fixed. One writes 

T(n, nu 2 ) = n n / — e -»(«-J»g«) = n ™ / dy e -nf(y) > = y _ logy ) ( B2 ) 

Ju 2 x J\ogu 2 

and under this form (|B2|) we can now perform a saddle-point calculation. The function f(y) has a single minimum 
at y = 1 and therefore we expect that the large n behavior of the expression in Eq. (|B2[) will depend on the sign of 
u — 1. For u < 1 the minimum of f(y) lies in the interval of integration in Eq. ()F32|) and one gets a result independent 
of u 

T{n + l,nu 2 ) = V2^n n+1/2 e- n (l + 0(7i- 1/2 )^) , t* < 1 . (B3) 

On the other hand, for u > 1 the minimum of f(y) in Eq. (|B2[) lies outside the interval of integration and one gets 
instead 

2 7/ 2 ' 1 / 1 II 2 \ 

r(n + 1, mi 2 ) = n"- x e-»" 1 - - . „ - - . + 0{ n - 2 ) , u>\. (B4) 



Using the asymptotic expansion for it < 1 (|B3[) together with Eq. (|B1[) one obtains straightforwardly that for u < 1 
fixed and large n 

Pn{Vnu) = - + 0(n- x l 2 ) , u < 1 , (B5) 

yielding, in the limit n oo the expression (|4^)l given in the text. Similarly, using Eq. (IB4[) together with Eq. (|B1[) . 
one gets for it > 1 fixed and large n 

p n (Vnu) = Oin- 1 ) , u > 1 , (B6) 
yielding, in the limit n — > oo the expression (j43[) given in the text. 
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APPENDIX C: COMPUTATION OF (N%([0, x])) c FOR KAC'S POLYNOMIALS AND d = 2 

1. Scaling form 

In this appendix, we give the details of the com put ation of ( N„ (T O, x])) which lead to formula (l80|) given in the text. 
We start from the general formula valid for all n [38| (see also 



f-X f-X f-X 

(NZ([0,x]))= dt Pn (t)+ dtx d£ 2 £„, 2 (£i,£ 2 ) , (CI) 
Jo Jo Jo 

where tC n , 2 (ti,t 2 ) is the 2-point correlation function of real roots of K n (x), given by 



OO pOQ 



lC n , 2 (t 1 ,t 2 ) = \ I I \yxy2\e-^ YQ "^ Y) d yi dy 2 , (C2) 

4?r 2 x /det A„(£i,£ 2 ) J-c 



00 j —00 



where Y = (2/1,2/2), A„(£i,£ 2 ) is the 4 x 4 covariance matrix of {K n {tx),K' n {tx),K n {t 2 ),K' n {t 2 )) and Sl„(£i,i 2 ) is 
the 2x2 matrix obtained by removing the first and the third rows and columns from A n (£i,£ 2 ) _1 . In the limit 
0<1 — £1 <C 1, < 1 — t 2 ^ 1 and n — > 00, keeping u\ — n(l — t\), u 2 = n{l — t 2 ) fixed one has 

1 „ — x 

(K n (h)K n (t 2 )) = ng{u x + u 2 ) , g(x) = 



x 



{K n (ti)K' n (t 2 )) = {K' n {t x )K n {t 2 )) = -n 2 g'( Ul +u 2 ) , 
(K' n (t 1 )K' n (t 2 ))=n 3 g"(u 1 +u 2 ). (C3) 

Using these relations (|C3|) . one obtains the matrix A„(£i, t 2 ) in this scaling limit as 

(ng(2u{) ~n 2 g'(2u 1 ) ng(ui+M 2 ) -n 2 g'(ui + u 2 )\ 

-n 2 g'{2u 1 ) n 3 g"{2 Ul ) -n 2 g'( Ul +u 2 ) n 3 g"(ui+u 2 )\ (r ,s 

ng( Ul +u 2 -n 2 ff '(u 1 + W2 ) ng(2u 2 ) -n 2 g'(2u 2 ) ' 

-n 2 g'(ui +u 2 ) n 3 g"(ui+u 2 ) -n 2 g'(2u 2 ) n 3 g" (2u 2 ) / 

from which one gets 

detA„(£i,£ 2 ) = n 8 det A(ui,u 2 ) , (C5) 

(g(2ui) -g'(2ui) g{ui+u 2 ) -g'{ui+u 2 )\ 

-g'{2u x ) g"(2u,) -g'{u 1+ u 2 ) 9 "{u } +u 2 ) 

,g(ui+u 2 -g'(ui+u 2 ) g{2u 2 ) -g'(2u 2 ) y > 

-g'(ui+u 2 ) g"{ Ul +u 2 ) -g'(2u 2 ) g"(2u 2 ) J 

From the structure of A„(£i,£ 2 ) in the scaling limit (|C4[) , one obtains the one of £l n (ti,t 2 ) in Eq. (|C2[) as 



il{ui,u 2 )- t}(.. „,\ n<„. „\) > (C8) 



n n (t 1 ,t 2 )=n- 3 n(u 1 ,u 2 ) , (C7) 

'A(ui,u 2 ) B(wi,u 2 ) N 
^B(u l7 u 2 ) C(u 1 ,u 2 ) / 

where Cl(ui, u 2 ), is obtained by removing the first and the third rows and columns from A(ui, u 2 ) _1 . The functions 
A(u\, u 2 ), B(ui, u 2 ), C(u\, u 2 ) have complicated expressions not given here which can be easily computed e.g. using 
Mathematica. Using the scaling forms (|C51 IC7[) one obtains /C„ j2 (£i, £ 2 ) in Eq. ()C2[) in the scaling limit as 

£n,2(£i, t 2 ) = n 2 iC(ui,u 2 ) (C9) 

(detA)~ 1/2 ( 5 \ B 

K= ) tttoT 1+ ; ArcSinf , 8= -== . (CIO) 

n 2 AC(l - ((5) 2 ) I y/T=W J 
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where we used the notation A = K(u\,u 2 ) and the expression [2(| (see also Ref. [381 ]') 



\yiy2\e-^ M+2By ^ +Cy ^d yi dy 2 



— oo J — oo 



: ArcSintf 



(Cll) 



Using this scaling form (|C9|) . together with the one for the density in Eq. (|28p. one obtains (iV„([0, x]) 2 ) c in the limit 
< (1 — x) <C 1) with n — > oo keeping y = n(l — x) fixed as 



(N n {[0,x}f) c ~ p K {u)du+ dm du 2 {t{u 1 ,u 2 ) ~ p K {u 1 )p K {u 2 )) . 

Jy Jy Jy 



Alternatively, one can write Eq. (|C12[) in the large n limit as 

(JV n ([0, x]f) - (N n ([0, l]) 2 ) = -i/_(n(l - x)) 



v-{y) = j dup K (u)+2 / dui / du 2 (IC(ui,u 2 ) - p K (u 1 )p K (u 2 )) 
Jo Jq Ju t 



(C12) 

(C13) 
(C14) 



2. Large argument behavior of v~(y) 

To compute V—(y) in the large y limit, one compute dv-(y) / dy from Eq. (|C14p 

dv-{y) 



dy 



= p K (y) + 2 du 2 (IC(y, u 2 ) - p K (y)p K (u 2 )) . 



The large y behavior of p K (y) has been computed previously (|30p . To extract the large y behavior of dv. 
thus needs to compute the behavior of K,(u\,U2) for u\,u 2 3> 1. We first analyse the behavior A(tii,U2) 
for ui,u 2 3> 1. This yields 



(C15) 

(y)/dy, one 
inEq. jC6} 



/ 2ui 



A(ux,u 2 ) 



i i 

Aui U1+U2 

1 1 1 

4wf 4-u? (M1+M2) 2 («l+'«2) 

1 1 1 1 

«l+«2 (tll+«2) a 2u 2 4«f 



\(U1+I 



U2) 2 (wi+112) 3 



1 

! » J 



(«1+«2)A 

/ 



1 

1 



from which one gets 



det A(iti, U2) 



ui - u 2 



256(uiU 2 ) 4 \Ul+U2 

From Eq. (|C16|) , one gets the behavior of the matrix f2 in Eq. (|C8[) as 



At \ a 3 I M i +"2 



B(ui,u 2 ) ~ 16 



mi - m 2 
u 2 it 2 ( U\+ u 2 



u\ + u 2 \ui — u 2 

4 



nt \ q 3 / "i + "2 
G(u 1 ,u 2 ) ~ 8«2 1 



which gives 5 = B/V AC as 



<5(wi,w 2 ) 



2yuiU2 

Ml + U 2 



(C16) 



(C17) 

(C18) 
(C19) 
(C20) 

(C21) 
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Finally, using Eq. (|C18IC21|) together with Eq. ^ in d = 2 one gets 

K.{u u u2)-p (ui)p (u 2 ) -, ■ + 2 ^ = 7 ■ ^ArcSm-^- . C22 

Therefore, from Eq. (|C15[) . one gets 

to-W- 1 ^ 1 Vo( y - 2 ), (C23) 



dy y \7r 2 7r 
where we have used the identity 



r^7^ A -Sin(^)=;-l. (C24) 
Ji sft{v + l) 2 \v + lj 2 



Finally, one gets from Eq. (|C23|) that 



APPENDIX D: COMPUTATION OF {N%([-x,x])) c FOR WEYL POLYNOMIALS 

1. Scaling form 

In this appendix we give the detail of the computation of (N n ([— x, x]) 2 ) which leads to formula (|86[) in the text. 
Here again we start from the general formula valid for all n, similarly to Eq. (|C1|) 

r-X r-X f-X 

(N^([-x,x\)) = / dt Pn (t)+ dh dt3W„, 2 (ti,ta) , 



(Dl) 

where W,i, 2 (ti, t 2 ) is the 2-point correlation function of the real roots of W n (x), given by formula (|C2[) where K n {x) 
is replaced by W n (:r): 



W«, 2 (ti,i2) 



47rVdet A„(ii,t 2 ) 



OO />OG -i 



(D2) 



oo « — OO 



where F = (j/i, y 2 ), A„(ti,t 2 ) is the 4x4 covariance matrix of (W n (ti), W^(ii), W n (t 2 ) 7 W' n (t 2 )) and fi n (ii,i 2 ) is the 
2x2 matrix obtained by removing the first and the third rows and columns from A„(ti, i 2 ) _1 - In the limit n — > oo, 
keeping t\ < y/ri, £ 2 < yfn fixed one has 

{W n {t{]W n {t 2 )) = e tlt2 , 
(W;(ii)W„(ta)> =t 2 e tlt2 , 
(W n (tiX(t a )) = tie t ^ , 

(^(iiXfe))-(l + tii 2 )e tlt2 , (D3) 



from which we obtain A„(i 1; i 2 ) = A(t 1; i 2 ) as 

A(ti,t 2 ) = 

The determinant is easily obtained as 







gtlt2 


tie* 1 ' 2 \ 




(l+t?)e*? 


t 2 e tlt2 


(l + tii 2 )e 41 ' 2 ) 
t 2 e' 2 


e tlt2 


i 2 e tlt2 


t 2 
e 2 




(1 +ht 2 )e tlt2 


■f 2 

i 2 e* 2 


(l + il)e' 2 / 



(D4) 



det A(tx,t a ) = e < tl+t2 ) 2 (4sinh 2 (^-^l!) - (*! - < 2 ) 4 



(D5) 
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From A(ti,t2) in Eq. (|D4[) . one obtains £l n (tt,t2) = ^(ii,^) for large n as 



with 



t a ) = ^TT(1 - e^^) 2 (l - (t t - £ 2 ) 2 )) , 

F(t!,* 2 ) = |^r(l - e-^-^'(l + fa - t 2 ) 2 )) . (D7) 

Finally, using these expressions (|D7[) together with the formula in Eq. (|C11[) . one obtains from Eq. (|E)2|) that 
W„, 2 (ti,t2)=H ; (ii-i2) with 

W(-) = -i ((1 - e 7 2)2 -f e " S2)l f 1 + - ^ (S) Anflin^ , (D8) 
^ 2 1 - e-^ 2 ^ y/l - (5{s)) 2 7 

with 



e - 2 / 2 (e- + S 2 -1) 

<5(s = — ; —2 2 _ S 2 ■ (D9) 

1 — e — s e 



Notice that VV(s) is the two point correlation function of real zeros of W n (x). Interestingly, this correlation function in 
Eq (ID8I) coincides (up to a multiplicative prefactor 7r~ 2 ) with the correlation of straightened zeros of B n (x) obtained 
in [20| | . Its asymptotic behaviors are given by 

fJ-L s ^ 

Finally, using the expression (|D1[) together with Eq. (|D2|) one has for < x < \fn, in the limit n — > 00 



(7V 2 ([-x,x])) c = ^(ar), 

v(x) = — + 2[ ds(2x- S )(W(s)-7r- 2 ) . (Dll) 



2. Large argument behavior of f(a;) 
To analyse the large x behavior of v{x) 1 one computes dv(x)/dx and gets immediately 



lim = - + 4 / ds(W(s) - it- 2 ) (D12) 



(i 



z^oo dx IT 

such that 

oc 2^000; , i>1 (D13) 
with i/qo = 0.181988..., which leads to the value of Ooc up to order 0(e 2 ) given in Eq. (f89|) . 
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APPENDIX E: COMPUTATION OF (N*([a, b])) c FOR BINOMIAL POLYNOMIALS 



In this appendix, we give the detailed calculation of (N^([a, b])) c which leads to the formula (|93|) given in the text, 
for m = 3. We start with the general formula (see for instance Ref. [20l|): 

(Nl{[a,b]))c = (N n [a,b}) + 3(N*[a,b]) c 



dti / dt 2 / dt 3 {B n , 3 (t x ,t 2 ,t 3 ) - 3B nt2 (tx,t 2 )p n (t 3 ) + Pn{h)p n {h)p n {h)) 



(El) 
(E2) 

where B n ^ m is the TO-point correlation function of real roots of B n (x). In Eq. (|E1|) . B nt2 (ti, t 2 ) is given by the formula 
(|C2p where K n (x) is replaced by B n (x) and B n ,3(ti, *2, is) is formally given by (see for instance Ref. [201 ]) 



Bn,z(ti,t 2 ,h) = / dyi / dy 2 / rfj/3|yiy22/3|-Dn,3(0, J/1,0,2/2, 0,2/35*1, *2, *3) , 



(E3) 



where D ni3 (xi, y\, x 2 , y 2 , x 3 , y 3 ) is the joint distribution density of (-B„(ii), B' n (t{), (B n (t 2 ), B' n (t 2 ), (B n (t 3 ), B' n (t 3 )). 
According to Eq. (|92|) . the two terms in (|E1[) have the desired form (|93l) for large n. To study the triple integral in 
Eq. (|E2j) in the large n limit, we will use the results obtained in Ref. [2fJ : 



lim 



B n , 2 (ti,t 2 ) 
Pn{tl)pn{h) 



ti =<iV[0, S( ]>- 



b 2 (si,s 2 ) = b 2 {si - s 2 ) 



where (N[Q, s,]) 1 is the inverse function of (iV[0, Sj]). Similarly 



lim 



S n , 3 (ti,i2> is) 
p n {ti)p n {t 2 )p n {t 3 ) 



U = {N[0, Sl ])- 



= b 3 {si,s 2 , s 3 ) = 6 3 (si - s 2 , s 2 - S3) 



(E4) 



(E5) 



where b 2 (u) and b 3 (u, v) are well defined functions, with a quite complicated expression not given here (see Ref. [201 ] 
for more detail). Given these results (|E4l IE5[) . it is natural to perform the change of variable Sj = (N n [0, U]} in Eq. 
p2|) . which yields 



dtl / dt 2 / eft 3 (B„ 5 3 (ti, t 2 ,t 3 ) - 30 n , 2 (*l, t 2 )p n (t 3 ) + Pn{tl)Pn(t2)Pn{t3)) 



r(K„[0,b}) 


/■<Ar„[o,6]) 


/•<W»[<M>]) 


/ dsi 


/ dS2 


/ ds3 (&3(si - s 2 , s 2 - s 3 ) - 3& 2 (si - s 2 ) + 2) 


/<Af„[0,a]) 


/<iV„[0,a]) 


/(JV„[0,o]) 



Performing the change of variables Si — » — (AT n [0, a]), one has 

dti \ dt 2 I dt 3 (Bn i3 (t 1 ,t 2 ,t 3 ) - 3B n , 2 (t 1 ,t 2 )p n (t 3 ) + p n (tl)pn(t2)pn(t 3 )) 



r{N n [a,b]) Miv„ia,oj; wjv»ia,oj; 

/ dsi ds 2 ds 3 (6 3 (si - s 2 , s 2 - s 3 ) - 3fc 2 (si - s 2 ) + 2) 

Jo Jo Jo 



(E6) 
(E7) 

(E8) 
(E9) 



Given that (N n [a,b}) oc ^/n in the large n limit for a,b > n 1 / 2 , one gets the multiple integral in Eq. (|E8[) in the 
large n limit as 



p(N n [a,b]) MJVnLO.oj; /■ 

/ dsi / ds 2 / 

Jo Jo Jo 



G?Sl 



<iV„[a,b]) /•<iV„[a,b]) 

ds 2 / ds 3 (b 3 (si ~ s 2 , s 2 - s 3 ) - 36 2 (si - s 2 ) + 2) 



(E10) 



with cr = 3 J"_ /_ dv (b 3 (u, v) — 3b 2 (v) + 2), provided this double integral over u, v is well defined (which we can 
only assume here given the complicated expression of 63(11, 1;)). 
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Finally, combining Eq. jET] E2 together with Eq. ((ET0|) one obtains that 

(N*[a,b]) c <x(3 3 (N n [a,b]) , (Ell) 

with 0,3 = 1 + 3a,2 + cr. Notice that the crucial point in the derivation of this relation is the fact that (N n [a, b]) oc ^fn 
for any fixed a,b > n -1 / 2 . One expects a similar mechanism to hold for higher values of m, yielding Eq. ([53]) in the 
text. 

(N™[a,b]) c <x$ m (N n [a,b]) . (E12) 
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